Measurement of phase fluctuations of Bose-Einstein condensates in an optical lattice 
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Even at zero temperature, there exist phase fluctuations associated with an array of Bose-Einstein 
condensates confined in a one-dimensional optical lattice. We demonstrate a method to measure the 
phase fluctuations based on the Fourier spectrum of the atomic density for a condensate released 
from the optical lattice. The phase variance is extracted from the relative intensities of different 
peaks in the Fourier spectrum. This method works even for high lattice strength where interference 
peaks disappear in the atomic density distribution. 
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I. INTRODUCTION 

For Bose-Einstein condensates in an optical lattice, the 
phase fluctuation is a significant quantity in the investi- 
gation of quantum phase transitions [ij. The transition 
from superfluid to Mott insulator is usually accompanied 
by significantly increased phase fiuctuations which can be 
manifested in the interference pattern of the condensate 
samples. The vanishing of the contrast of the interfer- 
ence fringes is widely regarded as a characteristic of the 
quantum phase transition. 

The simplest lattice configuration suitable for demon- 
strating phase fiuctuations is a one-dimensional (ID) 
standing-wave laser field loaded by Bose-Einstein con- 
densates (BEC). Such a ID optical lattice is actually 
a common tool to test quantum properties of the cold 
atoms in periodic potentials. It has been used to demon- 
strate phase coherence [1, Q , Bloch oscillations (3] , num- 
ber squeezed state 01, Josephson current Q, nonlinear 
self-trapping of matter waves 0, 0] , and so on. In theory, 
quantum fiuctuations in phase and atomic number are 
often illustrated by considering a BEC in a double-well 
potential '^-Tg']. Of course, experimental measurement 
of the phase fluctuation plays a key role in understand- 
ing the quantum process occurring in a lattice. 

In a pioneering work by Orzel et al. phase fluctu- 
ations of the subcondensates in a ID optical lattice were 
measured by using the interference pattern of the released 
condensates. The phase variance was extracted from the 
contrast of the observed interference peaks. However, 
at very high lattice depth, the typical interference peaks 
disappear completely due to the large phase fluctuations, 
and this method is thus not valid any more. In Ref. , 
it is also shown that, close to the Mott insulator, the van- 
ishing of the interference fringes makes it difficult to de- 
scribe the quantitative changes of the system controlled 



with further increased lattice depth. 

In this paper, we develop a method to measure the 
phase variance by employing the Fourier spectrum of the 
released atomic cloud. Particularly, a simple analytical 
expression is found to extract the phase fluctuations. Our 
method works in principle even when the visibility of the 
interference peaks is completely lost, as demonstrated in 
our experiment. It is expected that this method provides 
a unique tool to other phase transitions in cold atomic 
systems p^ . 

The paper is organized as follows. In Sec.|lll we present 
the theoretical model to extract the phase fluctuations 
from the Fourier spectrum of the atomic density for a 
condensate released from the optical lattice. In Sec. IIIIl 
the theoretical model is demonstrated in our experiment. 
The accuracy and validity of our method are discussed, 
and a comparison with the method in [H is also presented. 
Finally in Sec. IIVI we summarize our obtained results. 



II. THEORY 

We now consider a ID optical lattice in the dimension 
of the X axis, formed by a conventional standing-wave 
laser fleld. Its strength is usually measured in units of 
the recoil energy Ej. = h? /2m\'^ , where m is the atomic 
mass, h is the Planck constant, and A the optical wave- 
length. In the tight-binding limit, the condensate loaded 
to the lattice can be treated as a chain of disk-shaped sub- 
condensates equally spaced by the lattice period d = A/2. 
The total number of the lattice sites occupied by the con- 
densate is denoted by M . When suddenly released from 
the optical lattice, the condensate undergoes a free ex- 
pansion process. After a time of flight (TOF) of r, the 
wave function of the atomic cloud can be written as 
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where ^i{x,t) refers to the wave function of the subcon- 
densate initially centered at the ith lattice site, and \ai\ 
represents the probability for an atom roughly located at 
the Ith lattice site. The atomic density is then written as 

Lq 

M^l (1) 



n=l 



where Gn reads 



M 



1=1 



while Gn with n> 1 takes the following form: 
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a*iai+n^i{x, T)^i+n{x, t) + c.c. (2) 



Note that, Go is physically different from other G„ 
with n > 1. It contains no interference terms, and is just 
a direct sum of the atomic densities of all the subcon- 
densates. Therefore, it gives rise to a spatially smooth 
density profile. In contrast, G„ describes the interference 
between the subcondensates spaced by nd in the optical 
lattice. One characteristic associated with G„ must be 
mentioned. The integration of G„ is zero (J G„dx = 0) 
due to the orthogonality between different , which im- 
plies that Gn would give rise to interference structures 
rather than a smooth background in the density distri- 
bution. 



A. Fourier spectrum 

According to Eq. ([2]), the Fourier transform of Gn is 

Fn = -^y^aim+n I ^i{x,T)^l+n{x,T)e'''''<lx + C.C. 

v27r ^ J 

= -^Wn{k, t) V a^ai+ne'''^'^ + c.c. 
V2tt , 

(3) 

Here, 

Wn{k,T) = ^*{x,t)^{x - nd,T)e'-'"'dx (4) 



is independent of the lattice site I. The second line of 
Eq. ([3]) is obtained using the fact that all $/ are identi- 
cal wave functions except for their center positions. In 
the tight-binding limit the Wannier function ^{x,t = 0) 
can be well approximated by a Gaussian wave packet 
(ttct^)"-^/** exp(— a;^/2(T^), where a — yJhfnmT^ is the os- 
cillator length, m the atomic mass and uix/^t: the axial 
trapping frequency of the lattice wells. After the TOF, 



the expanded subcondensate has a much larger size than 
its initial wave packet (fir/m ^ tr^), then the wave func- 
tion of a single subcondenste can be written as [T9| 
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Substituting Eq. ([5]) into Eq. @ , one gets 



r) = W^„ / exp [ix (fc — nfci)] dx. 



(5) 



(6) 



where fci = 27r/Ai, Ai = 2'KhT/dm and 
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Here, Ai is a characteristic length, equal to the travel dis- 
tance of an atom with a velocity twice the single photon 
recoil velocity. The integration term yields a (5-function 
like peak at fc = nki , and the peak width is inversely pro- 
portional to the spatial size of the expanded wave pack- 
ets. Similar analysis to the conjugate part in yields 
an identical peak at the symmetric position, fc — — nfci, 
instead. Therefore, the Fourier transform of G„ shows 
a pair of peaks at fc = ±nfci (only one peak for Go). 
Apparently, the characteristic length Ai is actually the 
spatial period of the Fourier component corresponding 
to the peak of the n = 1 order. As such, coherence prop- 
erties associated with different spacings between subcon- 
densates in the optical lattice can be distinguished from 
one another just by inspecting the Fourier spectrum of 
the density distribution of the expanded atomic cloud. 
The whole power spectrum of the atomic density is sim- 
ply given by S{k) = E„ |F„|^. 

Note that, Ai is much larger than the initial condensate 
size, which means kid <^ 1, and hence e**^''' ~ 1. One 
may get from Eq. ^ the peak intensities in the power 
spectrum S{k): 



|F„(fc = nfci)r 



AYn 



(7) 



1=1 "/ ai+n 



where A = \wo{k — 0)| /27r and y„ 

In particular, Yq — 1, as required by the normalization 
condition. From the expression in Eq. ([5]) , one sees that 
the amplitude of Wnik = nfci), and hence A, is indepen- 
dent of n. Therefore, the relative intensity of peaks in 
the power spectrum depends only upon y„. 



B. Phase Fluctuations 

We now turn to consider the peak intensities in the 
power spectrum, from which the phase fluctuation can 
be deduced. In the optical lattice, the confined subcon- 
densates undergo phase fluctuations. The phase factor 
of each subcondensate is contained in the corresponding 
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coefficient ai, and the summation term y„ in Eq. ([7]) is 
then rewritten as 
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where (50;„ = <j)i+n — (f)i. Due to phase fluctuations of 
the subcondensates when trapped in the optical lattice, 
54>in takes random values with zero average. The last 
line of the above equation is obtained by assuming that 
|Q;;a/+„| changes slowly with I. As long as n <C M, the 
summation term |q;;Q!/+„| is constant for different n, 
i.e., |a;Q;;+„| ~ |a/|^ = 1. Yn is then simplified to: 



Yn = 



1 



M-n 



(8) 




FIG. 1: The solid line shows the numerically calculated Yn 
versus the variance of the random phase 54>in for a total lat- 
tice number M — 200. Each data point is an average of 1000 
individual runs of the computation. Shaded area corresponds 
to the standard error. The dashed line is an exponential curve 
in the form of e ^ . The inset is the same figure on a loga- 
rithm vertical scale, highlighting the discrepancy between the 
two curves. 



Apparently, for n > 0, y„ depends on the variance of 
Scfun, denoted by =< S(t>f^ >. As shown in Fig.[Tl 
Yn drops quickly with increasing ^. The randomness of 
S(j>in results in the fluctuations of F„ around its averaged 
values. In the region of larger phase variance, Y„ has 
larger fractional fluctuations. In computation, the total 
lattice site number is assumed to be M = 200. We find, 
however, the averaged value of Y„ is nearly unchanged 
for different values of M. 

Checking the summation e*"*"^'" in Eq. (|H]), one sees 
that it is nearly a real number as the imaginary terms are 



averaged to zero. By replacing S(f)in by its variance ^ and 
assuming ^ ^ 1, this summation term can be approxi- 
mated by (Af — n)e^ 2^ . We then have an approximated 
expression of Yn in an exponential form; 



Yn = e-« 



(9) 



This exponential curve is plotted in Fig.[l] as well. Sur- 
prisingly, it shows a good match to the numerical data 
(the solid curve in Fig.[T]) even for the region far beyond 
^ ^ 1. For example, the relative error is only 5% at 
^ = O.Stt. Of course, the tendency to larger relative er- 
rors can be clearly seen as ^ is increased. 

The phase fluctuation of two adjacent subcondensates 
can be measured by a phase variance — It is 

straightforward to prove that two subcondensates spaced 
by n times the lattice period correspond to a n times 
larger phase variance, i.e., = ncr^. One knows from 
Eq. dZl) and dH) that P„ = Ae^''-'''' . Taking natural loga- 
rithm for both sides of this formula, one obtains 



In Pn =hiA — na^ , 



(10) 



which clearly shows the linear relation between the loga- 
rithmic scale of the peak height and subcondensate spac- 
ing n. The slope of this linear curve is just the phase 
variance cr^. Therefore, the phase fluctuations of the con- 
densates confined in discrete lattice wells can be easily 
determined by the peak structures in the Fourier spec- 
trum of the atomic density distribution after the time of 
flight. We stress here that this method to measure 
does not need exact calibration of the peak intensities, 
because ct^ is unaffected by any identical scale factors 
applied to all P„. This clearly shows the convenience of 
Eq. (fTO)) in the measurement of cr^. 

There also exists an alternative method to measure the 
phase variance, where is determined by comparing the 
experimentally measured P„ with numerically calculated 
Yn- When ^ is very large, the analytical expression of F„ 
in Eq. Q is not a good approximation any more. In this 
case, the second method is more reliable than the for- 
mer one based on Eq. pUj) . We shall make a comparison 
between the two methods later. 



III. EXPERIMENT 



Our experiments are carried out by using a nearly pure 
^''Rb condensate in the hyperfine state \F = 2,Mp — 2) 
with typically 10^ atoms. The experimental setup was 
described elsewhere [l^. The ID optical lattice is formed 
by a retroreflected laser beam with a wavelength of 
A = 1064 nm. Its strength is calibrated using a method 
of Kapitza-Dirac scattering [17, 2Ql- The recoil energy 
is Er ~ h X 2.03 kHz. The lattice light is adiabatically 
applied to the cigar-shaped condensate along its axial 
direction during a time of 50 ms. After a holding time 
of 10 ms, the lattice light, as well as the magnetic trap, 
are suddenly switched off. Finally, an absorption image is 
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taken for the released atomic cloud after a 30 ms of TOF, 
by using a probe light directed perpendicular to the lat- 
tice beam. The experimental parameters correspond to a 
characteristic length Ai = 259 fim, which is much larger 
than the pixel size (9 //m) of our CCD camera. In princi- 
ple, the peaks up to the order of n = 14 can be resolved 
by the CCD camera. In order to obtain the statistics of 
the phase variance, the experiment was repeated at least 
eleven times for each lattice depth. 
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FIG. 2: Open circles show the power spectral intensity ob- 
tained from the density distribution of an expanded conden- 
sate released from an optical lattice with a depth of 34.6 E^. 
Four points at the top of three peaks are marked by crosses, 
and used to determine a phase variance of = 1.79(8) by 
a linear fit (dashed line). The peak close to fc/fci = 3 is 
caused by optical noises of the probe light. The inset shows 
the atomic density distribution of this atomic cloud along the 
direction of the lattice beam, showing clearly two side peaks 
due to the interference of the subcondensates. 

The power spectrum S{k) is obtained from an ab- 
sorption image as follows: Spatial frequency spectrum 
is calculated by Fourier-transforming the matrix of im- 
age. Taking the absolute square to obtain the 2D power 
spectral density, and S{k) is then obtained from the 2D 
matrix by making summation along the direction per- 
pendicular to the lattice beam. Thus, only the power 
spectrum along the lattice direction is remained in S{k). 
Finally, taking the natural logarithm of S{k), we can find 
the peaks and extract cr^ according to Eq. (|10p . 

Figure[2] displays a power spectrum corresponding to 
the absorption image of an atomic cloud released from 
the optical lattice with a depth of 34.6 Er- As predicted, 
S{k) consists of a series of peaks equally spaced by fci. 
The peaks with n — 0-2 are true signals of the cold 
atoms, whereas the peak close to k/ki = 3 is confirmed 
to be optical noise of the probe light itself. This noise 
peak appears occasionally in repeated experiments, even 
in the absence of the atomic cloud. 

A higher phase variance means a weaker phase corre- 
lation. In our experiment, subcondensate pairs with a 
spacing larger than 2c? can not yield a peak {n > 2) high 



enough to be visible in S{k). We thus inferred that the 
phase correlation of two subcondensates drops quickly 
with increased distance between them, which is a natu- 
ral consequence of the proportional relation of = na^. 




Position (pixels) 

FIG. 3: (Color online) (a) Power spectrum S{k) of an ex- 
panded condensate initially trapped at a depth of 68iJr. (b) 
A false-color absorption image of this released condensate, 
with a field of view of 0.55 mm x 1.63 mm. (c) Atomic den- 
sity in the x direction, which is obtained by integrating the 
pixels in height. 

To extract the phase variance cr^, we used the data 
points at the top of peaks in a power spectrum to per- 
form the linear fit in the form of Eq. ([TUl) . The standard 
error of a fitted slope (ct^) is usually small (< 5%). In 
contrast, the value of cr^ shows much larger fluctuations 
from shot to shot as the experiment is repeated under 
the same conditions. We thus take only the latter fluctu- 
ations into consideration when calculating the error bars 
of the phase variance. A linear fit to the three peaks in 
Fig. m yields a phase variance of cr^ ~ 1.79. At the cor- 
responding depth level of the optical lattice, the atomic 
density profile exhibits clear interference peaks along the 
lattice direction (see the inset in Fig. [5]). As the lattice 
depth is increased, the phase variance is increased, and 
the higher order peaks in S{k) with n > 1 are expected 
to become weaker accordingly. Figure[3] displays a typ- 
ical result in such a case. At the depth of 68 Er, only 
the peaks of rt = 0, 1 are observed. Higher order peaks 
are too weak to be identified. As seen from Fig. [3] (b) 
and (c) , the released condensate has lost the interference 
structures completely, in contrast to the side peaks in 
Fig. [2] In fact, despite the loss of interference peaks in 
atomic density distribution, the peaks in S{k) can still 
be observed for strong optical lattice, up to the highest 
level {120 Er) we have reached. 

The dependence of phase variance a upon the lattice 
depth is displayed in Fig.|4l The two sets of data (blue 
and red ones) were obtained using the two methods de- 
scribed in Sec. Ill Bl There exists a clear trend, where the 
deeper the optical lattice, the larger the phase variance. 
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The value of a grows from ~ 0.47r to ~ tt, indicating that 
the relative phase between adjacent lattice sites gains in- 
creased randomness, and tend to be a completely random 
phase. Below 50£'r, the blue data points are very close 
to the red ones. Beyond this level, however, the discrep- 
ancy between the two sets of data becomes larger with 
increased lattice depth. It can be simply understood by 
the fact that the exponential form of Yn (Eq. © ) is not 
a good approximation for large ^. Roughly speaking, 
the linear fit method can only be applied when ^ is less 
than O.Gtt, as is evident in Fig. 2] It is worthy to point 
out that, although the linear fit method is invalid in this 
case, the alternative method is still simple to extract the 
phase variance from the Fourier spectrum. 
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FIG. 4: (Color online) Measured phase variance a versus 
the lattice depth. For each data point, eleven independent 
runs were averaged. Blue data points were obtained based on 
Eq. (|10p . whereas, the red data points were obtained by com- 
paring the measured relative intensity of the peaks of n = 1 
with the numerical values of Y\. Black data points were ob- 
tained by fitting the atomic density profile using the method 
in [|. 

The accuracy of phase variance a is limited by at least 
two factors. First, as shown by the solid curve in Fig.[Tl 
Yn is less sensitive to ^ when ^ becomes large, particularly 
in a region around tt. Second, Yn itself is in fact a fluc- 
tuating parameter in principle, although we have taken 
its averaged value as the measured result. The fluctua- 
tion is evident from the fact that the relative intensity of 
each peak in S{k) varies from shot to shot. Therefore, 
in the deep lattice region, data points of a are always 
accompanied by large error bars. 

As a comparison, we have also extracted the phase 
variance by analyzing the visibility of the interference 
peaks of the released atomic cloud. This method was 
flrst demonstrated in ^], where the visibility is charac- 
terized by a quantity C, deflned as the ratio of the width 
of a single peak to separation between the peaks. The 
phase variance is determined by comparing the observed 
value of C with those obtained from simulated data sets. 
First, we fit the interference profile with three Gaussian 



peaks to get the value of C- Then, we calculate the in- 
terference profile by using a simple one-dimension model 
which is similar to that in our previous work [l^. Each 
subcondensate in a single lattice well is assigned a ran- 
dom phase. Those random phases are set in such a way 
that the phase difference between two adjacent lattice 
sites obeys a Gaussian distribution with a given variance 

. We convolve the calculated interference profile with 
a resolution function to account for the limited resolu- 
tion of our imaging system. To obtain the simulated C, 
we fit the convolved waveform with the same fit function 
applied to the experimental data. The simulation proce- 
dure is repeated many times to obtain an averaged value 
of the simulated C- 

The phase variances deduced from this fitting method 
are displayed in figure 4. It is obvious that, below a 
depth level of ~ 60i?r, the results of the fitting method 
agree well with our method. This further confirms the 
validity of our method. However, above 60 Er, the fitting 
method does not work due to the following reasons: To 
strictly follow the method demonstrated in 0], both the 
peak width and peak separation must be treated as fitting 
parameters. However, at high lattice-depth level, only 
one broad peak is left in the interference profile, as shown 
in Fig.|31^c). The fitting procedure cannot give reasonable 
peak positions. More specifically, the side peaks obtained 
from the fitting program deviate significantly from the 
positions where they should be located. In addition, the 
fitted width of the side peaks is not reliable. 



IV. SUMMARY AND DISCUSSION 

We have developed a method to measure the phase 
fluctuations of the subcondensates confined in a ID op- 
tical lattice. In our method, the Fourier spectra of the 
conventional absorption images of the released atomic 
clouds have been investigated. The phase variance be- 
tween adjacent lattice wells is deduced from the relative 
intensities of the peaks in a Fourier power spectrum. Our 
experimental measurements have displayed an increased 
phase variance as the lattice depth becomes larger, and 
also indicated that phase correlation of two lattice wells 
decreases quickly with the increased distance between 
them. Our method does not rely on the existence of 
interference peaks of the released condensates, and it 
works even for very large lattice depth. This method 
is a complementary to that demonstrated in Q , and will 
be a useful tool in analyzing phenomenons associated to 
phase fluctuations in optical lattice systems, particularly 
for the case close to the quantum phase transition. 

Our theoretical model is established in the tight- 
binding limit, where a condensate in the optical lat- 
tice can be regarded as a chain of subcondensates. We 
found that, for weak optical lattice, there is no multi- 
peak structure in S{k). The peaks start to appear when 
the lattice depth reaches a level of ~ IQE^.. Well re- 
solved peaks can be observed if the depth level is further 
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increased to ~ 20 Er- It sets a coarse boundary beyond 
which our model is apphcable. Below the level of 20 E^, 
one has to switch to the method in [5] to measure the 
phase variance. 

The optical lattice is not homogeneous due to the pres- 
ence of the magnetic trap which is used to support the 
atoms against the gravity. The harmonic confinement of 
the magnetic trap corresponds to a trapping frequency 
of 27r X 7.6 Hz, and it remains until the sudden release of 
the atomic cloud. For a total atomic number of 10^, the 
number of lattice sites that are populated is M ~ 200. 
Most atoms are distributed in the center region of the lat- 
tice where the tunneling rate J is nearly uniform. This 
assures the assumption that a is uniform over the opti- 
cal lattice. On the other hand, we did not see noticeable 
changes of a as the total atomic number is changed from 
4 — 15 X 10*. It is due to the fact that the tunneling 
rate is independent of the atom numbers in single lattice 



wells. 

In principle, our method can be extended to 2D and 
3D optical lattices by treating the power spectrum in 
each dimension separately. For a 2D optical lattice, one 
probe beam perpendicular to the lattice plane is enough. 
For a 3D lattice, however, an additional probe beam is 
required to detect the atomic density profile in the third 
dimension. 
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